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Abstract. The Shannon dimensionality of orbital angular momentum (OAM) 
entanglement produced in spontaneous parametric down-conversion can be probed 
by using multi-sector phase analysers pQ. We demonstrate a spatial light 
modulator-based implementation of these analysers, and use it to measure a 
Schmidt number of about 50. 
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1. Introduction 

The photon pairs generated by spontaneous parametric down-conversion (SPDC) in a 
non-linear crystal are correlated in their various properties, some of which can exhibit 
high-dimensional entanglement. The orbital angular momentum (OAM) of light is 
one such property, associated with phase structures of the form e lt ^, where £h is the 
OAM carried by each photon [2]. As £ is an integer and is theoretically unbounded, 
OAM offers a natural discrete space for exploring high-dimensional entanglement [3]. 
In SPDC with a Gaussian pump (with £ = 0), the spectrum of OAM correlations 
between the signal and idler photons is peaked at £ = 0, with tails towards high \£\ 
values. The width and shape of the spectrum, and therefore the number of OAM 
modes that constitute the two-photon entangled state, can be engineered directly by 
manipulating the structure of the beam pumping the crystal [4j [5j [6] or by tuning the 
phase-matching conditions in SPDC [7]. 

It is necessary to distinguish between the entanglement of the high-dimensional 
OAM state generated by the SPDC process, whose dimensionality is described by the 
Schmidt number K [8] , and that detected by the measurement stage of the system [1] , 
with Shannon dimensionality M [10] . M is in fact dependent upon both the generated 
down-conversion OAM spectrum and the finite detection capability of our analysers. 
The latter is expressed in terms of the number of modes D that an analyser has access 
to. By knowing the relation between these three quantities, we can infer K after 
measuring M with an analyser characterised by a known D. 

Determining the measured dimensionality M of tailored high-dimensional 
entangled states can be carried out by performing appropriate selective projective 
measurements [9]. One such measurement is based on pairs of multi-sector phase 
plates, placed in the two arms of a down-conversion system. Each plate has N 
azimuthal angular sectors, each of which introduces a 7r phase shift. The number and 
angular width of the sectors of the phase analysers placed in each of the signal and 
idler arms define the superposition of OAM eigenmodes in which the two-photon state 
produced by down-conversion is projected. By optimising the binary phase profile of 
the phase analysers, it is possible to maximise the Shannon dimensionality D of the 
measurement apparatus, for any number of sectors N of the two plates. 

By using angular phase analysers we infer the Schmidt number K, characterising 
the effective number of azimuthal entangled OAM modes. In contrast to previous 
works, which used micromachined phase plates [11] . we implement ./V-sector angular 
phase analysers using spatial light modulators (SLMs) [12] . Computer-controlled 
SLMs provide a fast, convenient and reliable way of producing holographic phase masks 
with arbitrary orientations and numbers of sectors, to be used in the measurement 
of the Shannon dimensionality of OAM entanglement. The use of multi-sector phase 
masks to probe high-dimensional states, as opposed to narrow single-sector analysers 
|13[ 114] , allows the measurements of tight angular correlations whilst maintaining high 
optical throughput. 

2. Theory 

2.1. Amplitude and phase masks 

The angular measurement of a light field can be achieved by employing an angular 
slit. The idea behind this approach is the angular analog of a linear slit: a linear 
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slit measures a field at one linear coordinate, with an uncertainty that depends on 
the width of the slit. An angular slit, on the other hand, measures a field at one 
angular coordinate with an uncertainty that depends on the angular width of the 
slit. As the angular position and orbital angular momentum form a pair of conjugate 
observables, the tighter the angular correlations, the larger the spread in the OAM 
observable [HJ [16] . If one seeks high measurement accuracy, the angular slit has to 
be very narrow and, in turn, this means that much of the light is blocked, yielding 
the problem of lower number of counts as the measurement uncertainty decreases. A 
practiced solution to this problem is to employ phase masks, instead of amplitude 
masks. 

The design of the phase mask consists of an angular step plate, which is 
characterised by a number of sectors N and by a set of angles that describe the 
position and width of each sector [TT]. Each alternate sector applies a it phase shift. 
The overall action of a phase plate of this kind is therefore to flip the phase of a light 
field in each sector, about the centre of the plate, and leaving the phase unchanged 
everywhere else. Note that the action of a phase plate does not affect the radial 
degree of freedom, as the design is radially invariant. This means that there is no 
coupling between different eigenmodes of the radial degree of freedom, which allows 
us to restrict ourselves to the azimuthal content of the measured state. This effect can 
be described in terms of OAM. A plane wave is turned into a superposition of different 
OAM eigenstates. The range of OAM eigenstates of which the superposition consists 
depends on the number of sectors and on their relative positions and widths Such 
effect is analogous to the effect of an amplitude mask, without the drawback of letting 
less and less light through as the angular uncertainty is decreased. 



2.2. Measurement of Hilbert space dimensionality 

We probe the two-photon state produced by down-conversion 

oo 

l*>= E «*W.®M>< (i) 

i=—oo 

described in the OAM basis, where \£) s and \£) i correspond to the states of the signal 
and idler photons respectively. By expressing the projection state associated with a 
phase-mask analyser oriented at an angle a as the superposition 

\A( a ))=Y J M\l)e Ua , (2) 
l 

the coincidence probability for a pair of analysers, oriented at a and /3 respectively, is 
given by 

P{a 1 P) = \{A{a) 1 B{P)\^)\ 2 , (3) 

where (A(a),B(j3)\ = (A(a)\ <g) (B(/3)\. The coefficients lt = |(^|A(0))| 2 = \\ t \ 2 
(with Y^ili — 1), defined by the profile of the iV-sector phase masks, determine 
the respective OAM spectrum. It is possible to design the arc sectors of each N- 
sector phase mask in a way that maximises the Shannon dimensionality D of the 
measured entanglement. For each N we used the optimal arrangement of sectors, 
based on results from TTJ. The maximum number of modes D that can in principle 
be measured by each of such optimal TV-sector masks can then be inferred from the 
theoretical distribution of eigenmodes as D = I/E^tI' as shown in [TT] . 
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Figure 1. Experimental setup implementing SLM- based coincidence detection 
with angular-sector phase masks (see text for details). The /3-barium borate 
(BBO) crystal is pumped by a 355 nm laser and imaged onto two spatial light 
modulators (SLMs). The outputs of the SLMs are then imaged onto single-mode 
fibres (SMFs). 



In contrast, the Shannon dimensionality D of each iV-sector mask used here was 
derived through a numerical model. The numerical model considered the distribution 
of the overlap between two identical A^-sector masks within a two-dimensional region 
with a Gaussian profile, as the orientation /3 of one of the masks was rotated with 
respect to the other, a. The maximum measurable dimensionality D for each optimal 
A^-sector mask is shown in fig. @] As already noted in [TT], it was found that, for N 
at least up to 14, D increases approximately linearly with N. Given the coincidence 
probability distribution P(a — (3) obtained from the numerical model, the Shannon 
dimensionality [1] can be directly calculated as 



where angle /3 is measured with respect to a. We implemented multi-sector masks with 
a number N of sectors between 1 and 16. The dimensionalities D obtained from the 
numerical model previously described and theoretical model, the latter obtained from 
the decomposition of each A^-sector mask into eigenmodes 7^ and the sum over the 
first thousand terms of each corresponding infinite series expansion |11) , were found 
to differ by less than 2.5%. 

3. Experimental results 

3.1. Experimental setup 

Implementing phase-mask analysers using computer-controlled spatial light modula- 
tors allows for quick and effective measurement of the Schmidt number K of the 
entangled state produced by SPDC. No optical elements need to be fabricated, physi- 
cally rotated or replaced when using a different multi-sector mask, as the measurement 
process simply involves displaying one of a set of different rotated Af-sector holograms 
on the SLMs and performing coincidence detection. 



D 



(4) 
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Figure 2. Typical best Gaussian fits of coincidence probability distributions 
P(a — (3) are shown for N = 5,8, 10 and cj> = 0. The inset shows the full 
— 180° < a — ft < 180° range used in the measurements, for the aforementioned 
values of N. Background subtraction was performed by assessing the experimental 
accidental coincidences. 

A 5 mm-thick /3-barium borate (BBO) non-linear crystal cut for type-I collinear 
SPDC acts as our source of entangled photon pairs. The crystal is pumped by a 
1 W UV laser to produce frequency-degenerate entangled photon pairs at 710 nm. 
The co-propagating signal and idler photons are separated by a non-polarizing beam 
splitter, and redirected to SLMs. The SLMs, onto which the crystal output face is 
imaged by a 2x telescope, are encoded with TV-sector phase holograms. The SLMs 
are then imaged onto single- mode fibres (SMFs), which couple the SLM output to 
single-photon photo-diode detectors (fig. [IJ, whose output is routed to coincidence- 
counting electronics. The coincidence counting has a timing window of 10 ns. Narrow- 
band, 2 nm interference filters are placed in front of the detectors to ensure that the 
frequency spread of the detected down-converted fields is small compared to the central 
frequencies. SLMs introduce great flexibility in our measurements, but this comes at 
the price of an overall lower detection efficiency, as the diffraction efficiency of SLMs 
is around 50 %. 

This detection configuration is insensitive to any overall phase factors. Therefore, 
while the conservation of OAM in the SPDC process would require placing mutually 
phase-conjugate iV-sector phase masks in the detection arms (i.e., the and tt phase- 
shifted sectors are inverted between the two masks), two identical phase masks can 
be used instead. These phase masks are self-conjugate in case of a tt phase shift. The 
finite pixel size of the SLMs places a restriction on the width of the sectors that can 
be displayed on the holograms. We show that we can implement optimal multi-sector 
phase masks with N = 1, ... 16. Suppressing the centres of the holograms, where the 
N angular sectors meet in a very limited spatial region of the SLM displays, did not 
turn out to be necessary. 

The phase-matching conditions of the down-conversion process for the BBO 
crystal were adjusted by slightly changing the orientation of the crystal with respect 
to the propagation direction of the pump beam |17j . This allowed us to increase the 
width of the orbital angular momentum spectrum, and thus decrease the width of the 
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Figure 3. Detected orbital angular momentum spectrum with projective 
measurements using traditional forked holograms. Shown is the coincidence 
probability P(£i,£i = —Is) for i s € (—30,30), for collinear (<f> = 0) and near- 
collinear (<j> = —2.3) phase-matching conditions. 

angular correlations distribution. The intensity profile of the down-conversion emission 
can be expressed as: / oc \E\ 2 oc sinc(0 + c£) 2 , where £ is the external emission angle in 
air, c = (|fc s | + \k.j\)L / {2n) 2 is a constant depending on the experimental parameters, 
and 4> = (\k p \ — \k s \ — |fcj|)L/2 (with L length of the crystal) determines the degree 
of non-collinearity of the process [18] . Measurements were performed for collinear 
(4> = 0) and near-collinear (ip = —2.3) phase-matching conditions. 

The coincidence probability distribution P(a — (3) was obtained by changing the 
orientation ft of the second phase mask over the range a ± 180°, where a is the 
orientation of the first. We found that, for the purposes of the experiment, a Gaussian 
distribution is an excellent empirical fit for the coincidence probability distributions 
(fig. [2]). The detected number of modes M, dependent on both the source and the 
detectors' properties, was obtained by substituting the Gaussian fit of the measured 
coincidences to P{a — (3) in eq. (j4|). 

3.2. Results and discussion 

The phase masks used have no radial structure. Any sensitivity of the detection 
apparatus to the radial quantum number p is therefore due to the spatial selectivity 
of the SMF coupling. Unlike full projective measurements of the down-conversion 
entangled state over a range of OAM eigenstates \£) (fig. EJ), a phase-mask 
determination of the dimensionality K of the source does not provide any direct 
information on the shape of the OAM spectrum. 

We consider the OAM spectrum generated by the source and the coincidence 
probability obtained from the numerical model, both fitted with Gaussian 
distributions, from which the dimensionalities K and D can respectively be obtained. 
We performed projective measurements of idler and signal over \£) and \—£) 
respectively to verify the validity of the assumption concerning the OAM spectrum 
(fig- IS])- The numerical model, from which the values of D for different N arc obtained, 
calculates the overlap in eq. ([3]) . Although the theoretical shape of the distribution of 
the coefficients jt, which characterise the action of a phase mask on a light field, 
has own distinctive features, its implementation on an SLM makes it possible to 
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Figure 4. Collinear phase-matching conditions. Phase mask dimensionality D 
from numerical model (blue points), measured M (red squares), and calculated 
dimensionality K (grey circles). We observe that D < K for N < 7; therefore, 
the calculated M saturates to D for any given N < 7. The solid grey line shows 
the best estimation for the number of modes K of the source, 35 ± 2. 
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Figure 5. Near- collinear phase-matching conditions. See caption of fig. [4] Given 
the wider spiral bandwidth in the near-collinear regime, K saturates to D for 
N < 9. The solid grey line shows the best estimation for the number of modes K 
of the source, 49 ± 2. 



approximate it with a Gaussian distribution. In fact, as the hologram representing 
the phase plate is rotated on the surface of the SLM, any imperfections (finite pixel 
size, surface roughness, slight unevenness of the phase shift, electrical fluctuations) 
influence the effect of the phase plate in a stochastic fashion. As this a very small 
overall effect, it does not change the dimensionality of the phase plate but rather 
smooths out the distribution of the coefficients. 

When the dimensionality K of the generated state is very different from 
the dimensionality D accessible to the detection system, the resulting measured 
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dimensionality, M, is given by the smaller of the two. For cases where K and D 
axe comparable, M can be approximated by noting that all three distributions are 
close to Gaussian in form and hence their contributing widths can be combined as: 

(5) 



which gives: 



M 2 ~ D 2 K 2 

M ~ DK (6) 
VD 2 + K 2 

from which the Schmidt number K can be derived: 

DM . (7) 
VD 2 - M 2 

Consequently, the source dimensionality K can be inferred from the theoretical 
coincidence probabilities and measured coincidence distribution. The calculated 
dimensionality for each TV-sector phase mask is shown in fig. U (collinear phase- 
matching) and [5] (near-collinear phase- matching) . We measured Schmidt numbers of 
35 ± 2 for collinear down-conversion, and 49 ± 2 in the near-collinear case. The results 
for source dimensionality obtained from the phase-mask analysis are compatible with 
the Schmidt number derived from projective spiral bandwidth measurement (fig. [3]), 
with the assumption of perfect single-mode detection. 

The mean visibility achieved in the experiment, defined here as the ratio between 
the mean baseline of the measured coincidence probability and the peak at a — f3 = 0, 
without background subtraction, is 90 % for collinear phase-matching, and 92 % for 
near-collinear. Systematic errors due to misalignment are found to be much larger 
than photon statistics uncertainties. 



4. Conclusions 



We have shown how multi-sector phase-mask analysers can be implemented using 
spatial light modulators, and used them to probe the effective number of modes in the 
high-dimensional bi-photon entangled state produced by parametric down-conversion. 
We used a set of several multi-sector analysers to infer the Schmidt number for different 
phase-matching conditions, and therefore, different widths of the OAM spectrum of 
the source. 
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